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Under the terms of this grant, studies were made of the structure 
of the plasmas pher a as. seen by radiosounder measurements aboard the 
Alovetti^li satellite. During the early phases of the work, magnetic 
tape data files were obtained from the NASA Ames Research Center to give 
a reasonably complete set of high latitude electron density profiles. 

Owing to internal data processing difficulties, Ames was able to provide 
these data only after considerable delay (& months) . 

During the period between the initial data request and data avail- 
ability, considerable effort was expended to develop models of ion flow 
in the topside ionosphere. Th^se models took both H and 0 into account 

and permitted various parameter studies to be made of the various factors 

+ 

which affect H escape in polar wind flows. The results of these studies 
were published as an initial paper in the 1972 "Critical Problems in 
Hague tospheric Physics" symposium held in Madrid, Spain* A copy of this 
paper is included here as Attachment I. 

When the NASA data became available, extensive computer programs 
were written to display the measured electron density profiles in ways 
useful to geophysical analysis. It was found that the expected mid-latitude 
trough was easily discernable in the nighttime ionosphere at locations ex- 
pected from similar observations of the plasmapause. In the dayside iono- 
sphere, however, it proved extremely difficult to find any trough-like 
phenomena in the electron density. Using the previously developed computer 
models, it was possible to study the region where the plasmapause appeared 
to be absent. ' It was found that over much of the dayside, large H + fluxes 
were computed well inside the plasmapause extending down to L-shells as 



as low as 2.5. This morning flow has been described in detail as a 
paper given as Attachment II which will soon appear in Planetary and 
Space Science. 

Further studies were undertaken to develop ideas concerning the 
relation between the location of the plasmapause as it is found in the 
equatorial plane and the location of the ionospheric trough. Briefly, 
it was found that a good correspondence can be expected only in the 
nightside during periods of moderate to disturbed magnetic activity. 

In the dayside the trough is at lower latitudes than the plasmasphere 
owing to the way H + must replenish magnetic field tubes of considerable 
volume. A paper discussing this effect has been submitted to the Journal 
of Geophysical Research and is included here as Attachment III. 



ATTACHMENT I 


BEHAVIOR OF THERMAL PLASMA IN THE 
MAGNETOSPHERE AND TOPSIDE IONOSPHERE 

Peter M. Banks 

Department of Applied Physics and Information Science 
University of California, San Diego 
La Jolla , California 92037 
U. S.A. 


Introduction 


Radio measurements of whistlers in the early 1950*8 [1] provided the first direct 
evidence for substantial amounts of thermal plasma (n ~ 50G/cni ) at altitudes of 
several earth radii. Using the charge exchange reaction of Dungey [2], Johnson [3] 
identified this plasma as H* of ionospheric origin. Further interest in the plasma of 
these very high regions was spurred by the spacecraft [4] and radio whistler [5] 
discoveries of the plasmapause, and the realization that the high density plasma 
envelope surrounding the earth was shaped in the form of a distorted torus [6]. 

Many of the early ideas about the plasmasphere tended to ignore obvious connection 
between the ionosphere and the thermal particles at large distances from the earth. 
Although it was often thought that the transition from the 0+ of the upper F 2 -region 
to the H + of the magnetosphere could be described in terms of quiescent equili- 
brium of a diffusive character, recent measurements [7] and theory [8] suggest the 
existence of far more dynamic flow states throughout much of the magnetosphere. It is 
this dynamic behavior which is the topic of the present review. Although the varieties 
of thermal plasma behavior are considerable, it now appears that many puzzling plasma 
features of the plasmasphere, the high latitude polar regions and the global topside 
ionosphere can be explained in a qualitative, and sometimes quantitative sense. 

To set the stage for this review. Figure 1 shows a recent model of the magneto- 
spheric environment. Of special interest for the thermal plasma behavior is the 
plasmapause, a sharp gradient in electron density aligned perpendicular to magnetic 
Isabella ranging from L = 4 in the morning sector to L ^ 7 in the afternoon. A 
schematic diagram of the equatorial projection of the plasmapause is given by Figure 
2, taken from the work of Chappell, et al, [9]. 

There appears to be little doubt that the principal factors involved in the 
formation of a plasmapause are the general cross tail magnetospheric electric field, 
the co-rotation induced terrestrial electric field, and the existence of a finite rate 
at which the ionosphere can provide thermal plasma to the magnetosphere. A simple 
model of the plasmapause can be constructed in the following way: The principal 

thermal ion of the topside ionosphere and magnetosphere is H 4 * which has its origin 
in the upper F 2 ~ region from the charge exchange reaction 

+ i + 

0 + H ==? H + 0 (1) 

v • 


with rates * 2.5 * 10 ^ T ** cm^/sec and k, « 2.2 * 10 ^ T ^ cm^/sec determined 
from recent experiments [10]7 1 

In equilibrium, (1) gives the H + /0 + density ratio as 

n(H + )/(n(0 + ) = 9/8 n^/nCOKyx^, (2) 
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When an excess of tt is present owing to inward ion transport. Equation (1) acts 
as a sink, converting into 0 + . When there is a deficit of H* resulting from an 
outflow of ions, (1) acts as a source and tends to reduce the 0* density ♦ These 
actions are closely related to the diffusive character of the upper ionosphere and 
virtually all production or loss of H* occurs in the altitude range 550-900 km. 

When unbalanced pressure gradients, gravity or other external forces are present, 
plasma flow is rapidly established along the magnetic field lines. In the absence of 
external electric fields transverse to B, the plasma flow is constrained to be parallel 
to 2 and the steady state density distribution along a field line reflects the 
bouyancy of If 1 " in 0 + with the lighter gas essentially floating .upon the heavier as a 
result of gravitationally induced electric fields parallel to B. 

As changes occur in the pressure and density of the 0 + of the F 2 ~region, 
variations propagate through the H* plasma and inward or outward flows may occur as 
the plasma system seeks a new equilibrium. For the most part, however, deviations in 
plasma density in the outer sections of the field tube are not large and a typical 
equatorial ion concentration is of the order 5 * 1G 2 c mT^ to 5 x 10* cra~-*. 

To progress further with the model, suppose that beyond a certain latitude (or 
L shell) a transverse electric field is present. If properly oriented, such a field 
will drive plasma away from the earth in a direction orthogonal to 2 and 2. If this 
field is sufficiently strong, the plasma will reach the magnetopause in a reasonably 
short period of time (6 to 12 hours) and will be permanently lost to the solar wind. 

The effect of this outward flow upon the ion density of the region exposed to 2 will 
depend upon the rate at which ionization can flow out of the ionosphere to replace 
the ions drifting across magnetic field lines* For H* this flow is greatly restricted 
by collisions with 0 + in the regions below 1000 km altitude and, as a consequence, the 
Ion density of the regions exposed to the magneto spheric convection electric field is 
very low (n e sl 10“ 2 ctn”^) [9j. Comparing the characteristic ion densities for the 
relatively stationary H* and that which is subject to E * S drift, it is seen that the 
plasmapause represents the approximate limit of penetration of the magneto spheric 
convection electric field. However, as discussed later, conditions other than 
2 x B drift can give rise to low H + densities in the transition zone of the topside 
ionosphere. In such cases, there Is no simple relationship between the plasmapause 
as seen far from the earth and the ionization discontinuities seen in the ionosphere. 

Making use of the foregoing ideas, it is possible to explain more clearly the 
meaning of Figure 2. In the tail of the magnetosphere, the convection electric field 
lies in the dawn to dusk direction and the thermal plasma flow is towards the earth, 
as indicated by the dotted plasma trajectories. Owing to the frozen field concept, 
there are two ways of viewing this motion. The first is to regard the plasma as being 
driven across fixed magnetic field lines. Plasma at any height above the equatorial 
plane is forced to move orthogonal to 2 and 2 which, in this case, implies not only a 
motion towards the sun, but also a velocity component towards the equatorial plane. 
Owing to the convergence of the magnetic field, there is a compression of the plasma 
which Is not uniform along a given field line and parallel redistribution flows are 
needed to maintain the plasma parallel force balance. Although it is not obvious in 
this formulation, the fact that the magnetic field lines are equipotentials insures 
that the E * B velocity of plasma across a given field fine varies in such a way that 
the plasma moves uniformly from one field line to the next. 

The second formulation of the E x B drift problem is based on this last condition. 
Flasma and magnetic flux tubes are considered to move together with the plasma con- 
tained in a given field tube undergoing compressional and redistributive effects as 
the field tube itself moves under the influence of the external electric field. Plasma 
contained in a given field tube can be lost or replenished only through the feet of 
the field tube; i.e„ t a dipolar field tube acts as a reservoir of plasma which is free 
to flow parallel to B in response to pressure gradients or other parallel components 
of external forces. 
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Returning to the diagram . shown in Figure 2, the inward flow of plasma from the 
1 represents both a compression of ambient plasma and a net equatorvard motion of 
p&asiaa in the ionosphere. As the magnetospheric flow passes the morningside of the 
easrth, no further changes in plasma density occur until the trajectories of the motion 
mxsve away from the earth towards the magnetopause, when a simultaneous rarefaction 
and poleward motion of the ionosphere occur. Although the night side compression and 
Wayside rarefactions are real and important effects, the overall ion density in this 
o*afcer convection region is very low owing to the very low plasma density in the tail 
o€ the magnetosphere . This, of course, is related to the relatively slow rate of 
iosxospheric plasma replenishment. 

Inside the plasmapause the magneto spheric convection electric field is mostly, 
fmt: not completely, absent. As a consequence, the thermal plasma of the plasmasphere 
tends to rotate with the earth, or at least with the atmospheric regions between 100 
t© 150 km where dynamo electric fields are generated. Owing to the influence of the 
convection electric field the idea of co-rotation is thought to he progressively 
pCBorer as one approaches the plasmapause 111]. An indication of the departure from 
co-rotation is given in Figure 2 by the black dots which portray hypothetical 
time marks In the equatorial movement of thermal plasma. If one traces the cycle of 
rotation between L » 3 and L = 4, it is found that there is a substantial slowing of 
plasma motion in the bulge region; i.e., a westward drift of ionization is present 
i/n the ionosphere as the field tubes are influenced by the convective electric field. 
A& lower latitudes such an effect is absent during magnetically undisturbed times 
«aa> d a more perfect co-rotation is achieved. 

f 

Diurnal changes in the F^-region produce flows of thermal plasma within the 
plasmasphere* In addition, the presence of an afternoon bulge in the equatorial 
pLlasma profile indicates that significant changes occur in the volume of the plasma- 
E^pheric magnetic flux tubes. In the sunlit sectors these volume changes are primarily 
expansions with a consequent cooling of the plasma and the establishment of an upward 
refilling plasma flew from the ionosphere. At night the opposite effect occurs. The 
TO>Ixme decreases, the plasma is heated, and H + is forced out both ends of the field 
rube into the nighttime p 2 -region. 

With the foregoing ideas providing a conceptual framework, it is possible to 
study a wide variety of experimental data to provide a quantitative verification. 

Each of the following sections is devoted to a particular aspect of thermal plasma 

behavior. 


Effects of H 


Flow 


A brief description of the chemical processes affecting the formation and loss of 
jtpd 6* in the topside ionosphere was given in the last section. Although the idea 
of chemical equilibrium at low altitudes matched by diffusive equilibrium above is a 
convenient concept, it bears little relationship to the actual flow conditions present 
as a consequence of ionization exchange between the ionospere and magnetosphere. 


The basic equations needed to describe the altitude distribution of 0 and H 
have been given elsewhere [12], Results of ion composition computations based on a 
realistic model of the topside ionosphere are shown in Figure 3. At low altitudes 
chemical equilibrium is assumed present. Owing to the increase in the ratio n(K)/n(0) 
with altitude, the H* density also rises rapidly. To introduce the possibility of H + 
flow, H* density profiles are shown for four values of the H density at 3000 km; 
l.e., 10-*, 5 * 10 , 8 x 10*, 9 x 10*, and 1.5 x 10 1 * ions/ezn. Each of these H + 
boundary densities corresponds to a different H + flux, the values for the given 
boundary densities Tanging from 1.35 x 10 8 ions/cm 2 /sec for the low H + density to 
- 2.3 x 10 ® Icois/cmVsec for the high H + density. Since a negative flux corresponds to 
plasma inflow, the rightmost profile in Figure 3 is typical of a field tube which 
Supplies the ionosphere with H+ and, throuch charge exchange, 0 + . In a similar manner, 
the lowest IT*" densities are reached when flows out of the ionosphere. The condition 
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of diffusive equilibrium is not shown, but lies beweeen the two closely spaced 
profiles. The fact that the H + flux changes from 2.5 x 10' ions/cm 2 /sec to -2.1 * 10® 
ions /cm 2 /sec for a 12% change in the H + boundary density gives an indication of the 
.sensitivity of the ionization flux to the high altitude H + density. 

With regard to the 0 + density, a range of values is obtained which depends upon 
the H+ profile. The most rapid decrease of 0 + density with altitude in Figure 3 Is 
associated with the highest density H + profile. Likewise, the least rapid decrease 
-of 0 + density with altitude is found for the low density profile. Such a behavior 
relates directly to the magnitude of the gravitational induced charge separation 
electric field acting parallel to $. When H + is the dominant ion, this electric field 
Is small and 0*, as a minor constituent, receives little support so its density de- 
creases rapidly with altitude. When H + is a minor constituent, a large charge 
separation field is present and the Q + density falls off slowly with altitude. 

For many purposes the interchange of ionization between the ionosphere and mag- 
netosphere is best described in terms of the reservoir model given in the last 
section. Within the plasmasphere , field tubes daily expand and contract under the 
Influence of electric fields, thereby inducing continual ionisation flow. The mag- 
nitude of the If 1 " flux associated with changes in the reservoir (3000 km) density is 
illustrated in Figure 4 using the same ionospheric model shown in Figure 3. In this 
model the H + density at 3000 km is greater than 9 * 1 q 3 ions/cnr plasma flows into 
£he ionosphere, while for smaller densities an outflow occurs. The extreme sensi- 
tivity of the inflow flux to reservoir density is apparent. 

f 

A different behavior is found as the boundary denstiy is lowered since the 
fl + outflow rapidly approaches a limiting flux such that large decreases in density 
produce no change in flux. Such a situation is not related to the subsonic or super- 
sonic nature of the flow (the density profiles of Figure 3 correspond to subsonic 
flow) but depends upon the rate of diffusion through 0 + (the barrier effect) and 
the rate of H + production via reaction (1) [13] . 

Since changes in the H + reservoir density affect the E + flux, it is to be ex- 
pected that changes in the ionospheric Fj-region 0 + density at 500 km will produce 
similar effects. To illustrate this. Figure 5 shows the variation in flux when 
the 3000 km H + density is held fixed at 10 3 ions/cm 3 while the 500 km 0 + density is 
changed. Diffusive equilibrium for the coupled H^/0^ system occurs when the 0 den- 
sity is about 5 * 10 4 ions/cm 3 . If the 0 + density is larger than this value, there 
will be an outflow of while a smaller value leads to an inward flow. Results such 
as these are clearly important to the question of the nighttime maintenance of the 
f^-layer [14]. 

- An important aeronomical aspect of thermal plasma flow relates to the density 
of the neutral hydrogen needed to produce any H + flow. The dependence of the If*' flux 
on this parameter is shown in Figure 6 where three neutral hydrogen densities have 
been used with the previous 0 + boundary value of 1.38 * 10 3 ions/cm 3 at 500 km. As 
hydrogen becomes less abundant, the diffusive equilibrium H + boundary density 
diminishes. As predicted by theory, the limiting flux is found to be directly 

proportional to the neutral hydrogen density 113] • Since 7 x 10 4 atoms/cm 3 is not 
atypical of the hydrogen concentration for moderate solar activity, it appears that 
the rate of plasma flow will have a significant variation between solar minimum and 
solar maximum conditions. An example of the short term changes in neutral hydrogen 
density has recently been given by Mayr and Brin ton [15]* 

The altitude distribution of H* depends strongly upon its flow speed relative^to 
0 + and the relative H + /0 + density ratio. Examination of Figure 3 shows that once H* is 
the dominant ion there is no apparent difference between the inward or outward flowing 
solutions. Such a result follows from the basic equation of motion for R + and helps 
to explain why density distributions which seem to be representative of diffusive 
equilibrium can actually correspond to substantial ionization transport. 
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To show this, the steady state H + momentum equation can be written as 
du 1 dp 1 dp 

u — + — + ^ + g| f = - v u (3) 

ds n m ds n m ds 
e 

_v — V 

where u is the transport parallel to B, s Is a coordinate along B, n and n e are the H 
and electron densities, m is the HT^ mass, p ** n k and p e = n 0 k T fi are the and 
electron gas 'pressures, is the parallel component of gravity, v is the h + /0 + trans- 
port collision frequency, and the 0 + is assumed stationary. 

Consider the high altitude case where H + is the dominant ion; i.e., n « n fi . 
Equation (3) gives 


du 1 dp 

u — — + “ ■ * + g = — v u (4) 

ds n m ds 
e 


where p 

P 


P e + P* 


When H is the dominant ion far above the F 2 -layer, the term v u of (4) is nor- 
mally very small. If the flow is such that the term u du/ds is small, (4) can be 
written as 


1 <JP S 

^ IT + 8 „ - °’ 


<5) 


precisely the distribution attributed to diffusive equilibrium, even though u itself 
may be large and of either sign. 


Another feature of (5) involves the importance of thermal gradients in determining 
the actual density profile- Letting + T^, (5) can be written as 


1 dn mg 1 dT 

« +131 + 0 

n ds kT T ds 

e P P 


( 6 ) 


so that when dT /ds ^ mg M /k thermal gradients, rather than gravity will determine the 
altitude variation of the density. For mid -latitudes, thermal gradients greater 
than 0.5*K km will significantly affect the H* altitude distribution, a positive 
gradient reducing the density more rapidly with altitude than would gravity acting 
alone. 


Thermal Plasma in the Plasmasphere 

The diurnal exchange of plasma between the topside ionosphere and plasmasphere has 
been indirectly measured by incoherent scatter radars at Millstone Rill, Massachusetts 
and Arecibo, Puerto Rico [16]- The Millstone Hill results show large inward and out- 
ward flows on the order 10$ ions/cm 2 /sec during magnetically quiet perids; i.e-, when 
the plasmasphere is probably in a steady state condition. 

The conditions establishing the diurnal flow are found both within the ionosphere 
and the magnetosphere. Even if there were a strict co-rotation of the plasmasphere, the 
changes in F 2 -region density and pressure are large enough to drive large H* flows in the 
high altitude regions. For example, between 0600 LT and 1200 LT it is not uncommon for the 
0"*" density at 500 km to increase by a factor of five or more. Referring to Figure 5, 
such an increase leads immediately to the flux limited flow of H* which, owing to the 
large volume of mid-latitude magnetic flux tubes, cannot greatly affect the 3000 km 
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reservoir boundary density for a long period of time* (If 10 Ions /cm /sec flow Into 
a volume of cm , a period of 10^ sec is needed to increase the If*’ density from 

10^ ions/cm^ to 2 * 10^ ions/cm^.) This implies that for this simple model the plasma 
density of the plasmas phere is determined by the relatively long filling time and 
should not change greatly from day to day since any surplus reservoir density will be 
lost very rapidly (there is no limit to the rate of inward flow) * 


Although such a simple model is useful conceptually, it does not adequately 
describe the actual variations known to exist within the plasmasphere. To begin with, 
the Fj-region undergoes continual variations which are related to the presence of 
thermospheric winds- Since these changes are of a somewhat irregular nature, the 
ionization of the plasmasphere must undergo a continual agitation with the presence 
of compression and rarefaction waves* 


Perhaps a more serious variation is related to transverse electric fields having 
their origin in the mid-latitude dynamo system [17] and the partial penetration of the 
magnetospheric convection electric field into the outer plasmasphere [18]. The dynamo 
electric field is established in the lower thermosphere by the flow of neutral winds 
driven by solar heating of the neutral atmosphere. Theoretical studies have shown 
that^the electric fields established by the dynamo system result In certain regular 
£ * B drifts which prevent a strict co-rotation of the high altitude plasma- Typical 
values of are of the order 1 to 5 m volt/m can give ionospheric plasma drifts as 
large as 100 m/$ec. Such drifts, if directed towards the magnetic poles, imply 
changes in the volume of the magnetospheric magnetic flux tubes and a redistribution 
of enclosed plasma; I.e., flow between the ionosphere and the plasmasphere. 

A more Important effect of the same type arises from the partial penetration of 
the convection electric field discussed in the first section. Preferring to Figure 2, 
magnetic flux tubes between L = 3 and L — 4 tend to be slower than co-rotatian in the 
afternoon bulge region and faster in the morning sector. In addition, field tubes in 
this same region have simultaneous changes in L value which can greatly alter the 
total volume* For example, a field tube at L = 4 in the morning sector may expand to 
L . » 6 in the bulge region- (A schematic sketch of this effect is shown in Figure 7.) 
Since volume * there is a five fold increase in volume, a similar decrease in 
icn density, and, if therelation PV Y holds, a factor of 15 decrease in the total 
plasma pressure. Such a change in the plasmasphere pressure must lead to a replenish- 
ment flow of H* with a consequent decrease in density in the topside transition 
region; i.e., the flow conditions of Figure 3 apply. 

In the nighttime an analogous compression of high altitude plasma must occur with 
there being a large increase in plasma density and pressure as the field tubes drift 
from i * 6 to L ■ In fact, such a compression of an already dense plasma can lead 
to a sharper density gradient at the plasmapause than might otherwise be expected. 

Although it is convenient to assume that motions of field tubes are accompanied 
by a very rapid redistribution of the thermal plasma along the field tube, such a 
process may be sufficiently slow to permit sizeable parallel pressure gradients to 
exist. In the nightside compression, for example, there Is a differential effect 
which tends to compress the plasma more strongly in the equatorial region than at other 
points along the same field tube. This creates a positive pressure gradient which 
drives an initial flow towards the earth at speeds of 1 to 10 km/ sec. 


It is clear from the foregoing discussion that drift motions across L shells will 
drive plasma ioto or out of the plasmasphere. While the inward flow rate has no 
theoretical limit, the outflow does and in this situation the 0 + to H + transition 
moves upwards to great heights (see Figure 3). 


Effects of Magnetic Disturbances 

Observations of the location of the plasmapause show a strong dependence upon 
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magnetic activity [19]- This variation is clearly shown in Figure 8, taken from the 
work of Chappell et al. I9J and based on OGO-y thermal ion measurements in the night- 
time regions. In essence, increases in the level of magnetic activity move the 
plasmapause to lower latitudes through an intensification of the convection electric 
field. The thermal ionization in the newly exposed field tubes is rapidly swept away 
and an upward flow of H* from the ionosphere immediately begins; i.e. , the polar wind 
Is established to the new- plasmapause* 

As magnetic activity declines, the convective electric field weakens such that 
plasm is no longer swept away to the magnetopause over the region previously exposed. 
In other terms, the field tubes once again co-rotate with the earth and are subject to 
the diurnal ebb and flow effects. However, since virtually all plasma in the field 
tube was lost during the disturbance, a long period of refilling mist follow before a 
steady state is reached. In the initial refilling period the thermal plasma density 
is sufficiently low that the upward H + flux Is supersonic I20J . This implies the 
formation of a shock travelling down the field tube and an eventual transition to 
subsonic flow* The upward flux of ionization is unaffected by these matters, however, 
and continues at its flux limited value until nearly steady state conditions are 
reached. 

Observations of the refilling process have been made 121] using the characteris- 
tics of radio whistlers. Owing to the relatively slow rate of refilling, it appears 
that 5 to 8 days are required to replenish the emptied mid-latitude field tubes. 

Since magnetic disturbances tend to occur more frequently than this, the plasmasphere 
recovery from one storm may not be completed before another cycle starts. This would 
tend to keep the exposed regions in continual flow and reduce the plasma density below 
equilibrium values. 

Another aspect of magnetic storms Is the way in which the high altitude plasma 
of the old plasmasphere is lost. On the day side of the earth a stronger convective 
field results in a drift of plasma towards the magnetopause (towards high latitudes 
in the ionosphere). This Is equivalent to a reduction in plasma pressure and the 
Ionosphere responds with an upward H* replenishment flow. On the night side there is 
a tendency to compress the plasmasphere towards the earth; i.e., plasma is moved in- 
wards to lover L values and lost through parallel flow through the feet of field tubes 
to the upper region. This flow will enhance the F 2 ~layer density and move the peak 
to a higher altitude. Again, it must be realized that there is a clear relationship 
between cross 1 shell movements of ionospheric plasma and an associated flow of 
plasma along the lines of magnetic force. 


Thermal Plasma Behavior at High Latitudes 


Several recent reviews have been given of the conditions affecting the topside 
thermal plasma in those regions of the earth lying outside the plasmasphere {22], 

Owing to the presence of the convection electric field plasma is continually swept 
away from the earth in the direction of the front of the magnetopause (see Figure 2). 
This sweeping effect is responsible for electron and ion pressure gradients directed 
along the magnetic field lines. This, in turn, establishes an upward flow of H* whose 


source, given by (1), lies in the upper I^-region. Because the sweeping action of the 
E x b drift is rapid, the plasma densities within the magnetosphere are low and the 
ionosphere supplies fi at the maximum rate possible with fluxes in the range 7 * 10 7 
to 2 « 10® ions 


sec 


Typical ion density profiles computed for a summer polar cap are shown in Figure 
9. The corresponding IT velocity, given In Figure 10, indicates that the flow 
eventually becomes supersonic. Although the term polar wind was originally adopted 
[23] to describe this supersonic flow for the regions where the geomagnetic field 
lines connect to the magneto sheath, theory and experiment show that virtually all 
regions outside the plasmapause have large escape flows of As shown previously. 
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the magnitude of the flux is virtually independent of the supersonic condition. 

Hence, the term polar wind can be used to describe the general high latitude plasma 
escape process, even though the. flow itself may not be supersonic at a given altitude. 

To illustrate the geographical extent of the polar wind, data obtained from 
Explorer 31 measurements 124] are shown in figures 11 and 12 to represent summer and 
winter conditions. Each solid line represents an observation of greater than 1 km/sec 
H* flow, while the dotted lines indicate periods when no flow was seen. The fluxes 
were found to be of the order 5 * 10^ ions/cm^/sec in the summer and 10® ions/cm^/sec 
in the dark winter polar regions. The difference probably arises from a change in 
neutral hydrogen density associated with different temperatures of the thermosphere. 

Independent evidence for the polar wind flow outside the plasmapause has recently 
been found in high altitude Alouette II electron density profiles. To illustrate, 
profiles of electron density as a function of invariant latitude are shown for five 
altitudes in Figure 13. This nighttime pass shows the presence of a sharp plasmapause 
near 60°A. Inside the plasmasphere the electron density contours are closely spaced, 
indicating a large electron density scale height and the presence of H + in substantial 
quantity. Poleward of the plasmapause a more erratic behavior is seen with the high 
altitude variations between 67* and 74°A representing auroral zone heating. Even in 
this region, however, the electron density contours give plasma scale heights of the 
order of 400 to 600 km; i.e. , typical of alone. Detailed analysis of records such 

as these has shown that the presence of 0 to 2500 km can be explained only if the 
neutral hydrogen density is very low; i.e., much less than 10^ cm“^, or if there is a 
polar wind flow of with a flux in the range 7 * 10^ to 2 x 10® ions/cm^/sec. Since 
measurements 115] rule out the possibility of the very low neutral hydrogen densities, 
it must be concluded that the observed behavior is a consequence of H* flow. Another 
example of this behavior at high latitudes is shown in Figure 14. 

The schematic trajectory of plasma drift at high latitudes can be obtained from 
Figure 2 through projection of the dipolar field lines into the ionospheric regions. 
Figure 15 shows a typical drift pattern at high latitudes including both the plasma- 
pause and the extension of the polar cusp. 

Although it is convenient to think of thermal plasma convection in terms of a 
steady motion, numerous experiments have found that convection is strongest during 
polar substorms 124]. While the convection patterns seem to be consistent with the 
original theory of magneto spheric convection, the plasma drift is at times greatly 
enhanced. To illustrate, Figure 16 shows the results of a 24 hour measurement of 
plasma drift velocities at Chatanika, Alaska (L ■* 5.6) obtained by incoherent scatter 
radar. The data, plotted with the geomagnetic pole at the center of the diagram, 
show that in the late afternoon and early evening the ^ x B flow velocities are 
small. With the onset of a substorm near 2000 LT a rapid southwesterly drift begins. 
Near midnight there is a reversal and the flow is almost due nast {magnetic). During 
the rest of the day the flow is generally directed back over the polar cap. Comparison 
of these results with the schematic flow trajectories of Figure 15 shows many similar- 
ities. Thus, at least in these disturbed times it can be concluded that the large 
scale convective field needed to sweep high altitude plasma to the magnetopause is 
present. Similar measurements made during periods of low substorm activity have shown 
that the nighttime enhancement of drift velocities is largely absent and the observed 
drift velocities remain less than 150 m/sec throughout the day and night. Although it 
might be thought that the slowness of convective drift would lead to a rapid filling 
of the polar magnetosphere with thermal H*, the time scale for this to happen is of 
the order of several days, substantially greater than the time between substorms. 
Consequently, the refilling process is never able to create a substantial H* concen- 
tration in the topside ionosphere outside the plasmasphere. 

As illustrated by Figure 1, two types of field lines lead from the high latitude 
ionosphere: those connected to the -magneto sheath and those forming closed dipole-like 

loops. When the supersonic polar wind is convected onto a closed field line a shock 
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front is formed which propagates down the field tube towards the earth. Calculations of 
the time needed for this shock to reach the ionosphere and convert the flow to sub- 
sonic speeds indicate periods of 10 to 20. hours. To illustrate, a plot of plasma 
density as a function of distance along a convec ting magnetic field tube is shown in 
Figure 17. Although shock, transition -moves towards the earth rapidly at first, it 
eyentually slows so that the actual transition to subsonic flow takes an extended 
period of time. Since the upward plasma flux is unaffected by this transition, it 
appears that the polar wind flow -must continue as a general feature for both disturbed 
and undisturbed conditions. [25] 


Future Studies 


Our knowledge of thermal plasma behavior in the magnetosphere is based on a 
relatively small number of experiments . Many important phenomena thought to exist; 
i.e, , reconnection shocks, flow of plasma in field tubes, and plasma compression and 
rarefraction have yet to be confirmed. Nevertheless, it is clear that the thermal 
plasma of the topside ionosphere is not quiescent and must undergo a wide variety of 
motions consistent with changes in the underlying ionosphere and the surrounding mag- 
netosphere. 

Within the plasmasphere the most important study relates to the validity of the 
plasma reservoir concept; i.e., are cross-L drifts of plasma accompanied by inward or 
outward (parallel to 5) flows? 

Although many studies assume a diffusive type density distribution of plasma along 
lines of magnetic force, proper measurements have never been made and the actual 
density profile may show pressure gradients associated with shock fronts or compression 
and rarefraction effects. 

Outside the plasmasphere virtually no magneto spheric thermal plasma measurements 
have been made to determine parallel flow speeds, temperatures, or ion flux. 

In many respects theory leads experiment in the investigation of thermal plasma 
behavior. Before further significant progress can be made it will be necessary to 
obtain a better idea of the actual physical conditions present in the various regions 
of the magnetosphere. 
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Captions 


Figure 


Figure 


Figure 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


Figure 8 


Model of the Earth’s Plasma and Magnetic Field Environment. The dark arrows 
indicate polar wind ion flow outside the plasmasphere. Similar flows some-- 
times present inside this region are not shown. 


Plasma Drift Motions Viewed iu the Magnetic Equatorial Plane. Outside the 
plasmapause the convection pattern lead to a loss of thermal plasma at the 
magnetopause. Inside the plasmasphere co-rotation effects are most important- 
The relative speed of this plasma is indicated by the black dots. Because 
the magnetic field is of a dipole nature, there is a considerable change in 
field tube volume implied- in all cross - L motions. 


Topside Ion Composition witjji H Flow Effects. The model is based on a 1000°K 
model thermosphere and an 0 density of 1.38 x 10-* ions/cnr at 500 km. Chemical 
equilibrium is assumed at 500 km and ti flow is obtained by varying the H* 
density at 3000 km. The numbers above each profile give the H + flux at 3000 km. 
A positive flux corresponds to upward flow. A neutral hydrogen density of 
1.92 x 10*/cnr w 3s chosen at 500 km, 

+ + 

. H flux for Different H 3000 km Boundary Densities. The fluxes of Figure 3 
are plotted as a function of the boundary density. Diffusive equilibrium 
corresponds to zero flux. Large H + inflow follows for relatively small 
increases of density. The H* outflow has a limiting value. 

4 > - 4 * 

. Variation of H Flux with 0 Boundary Density. The H density at 3000 km 
is fixed at 10^ ions /cm and the 500 km 0 + density is varied. Diffusive 
equilibrium for H occurs at n(0~^) — 5.5 * 10^ ions/cnr^. Smaller values 
lead to large inflow while the outflow increases gradually as larger values 
are chosen. 


. Variation of H* Flux with Neutral Hydrogen Density. The model given in 
Figure 3 was used with different neutral hydrogen densities. Both the 
limiting flux and the equilibrium H + density are affected by these changes. 


Schematic Variation of Magnetic Flux Tube Volume. As field tubes move 
across L-shells their volume changes as Volume * L^. Plasma redistri- 
butions are brought about by pressure gradients which move plasma along 
the field tube. 


. Location of the Plasmapause. Experimental OGO-V results showing the 
location of the plasmapause for different levels of magnetic activity. 
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Figure 9. High Latitude Ion Composition Profiles. Theoretical models were used 

to compute the ion composition for the polar regions. The flow becomes 
supersonic above 4100 km. The large H + flux is a consequence of a 
slightly high neutral hydrogen density and a large 0* density (see Figure 
5). 

Figure 10. Profiles of High Latitude H + Velocity. The profiles are given for three 
model neutral atmospheres and identical F^-region photoionization rates. 
The crosses indicate the point of transition to supersonic flow. 

Figure 11. Experimental Observation of Summer H + Flow from Explorer 31. The 

arrows indicate portions of the satellite path when data was obtained. 
Dotted arrows are plotted for times when no flow was seen (< 1 km/sec). 
The extent of the solid arrows support the idea of plasma flow in 
all areas outside the plasmasphere. 

Figure 12. Experimental Observation of Winter H Flow from Explorer 31. Owing to 
the very low plasma densities relatively fewer observations were 
possible. Nevertheless, the ion flow was found in the expected regions 
outside the plasmasphere. 

Figure 13. Observations of High Latitude Electron Densities from Alouette II. The 
topside radio sounder provided altitude profiles of electron density. 
These summertime plots show contours of electron density as a function 
of invariant latitude in the midnight sector. The plasmapause is 
seen in the topside near 6Q°A while auroral zone heating centers at 
70°A. The contours indicate only 0 + outside the plasmapause; i.e.» 
a polar wind flow is present. 

Figure 14. Observations of Summer High Latitude Electron Densities from the 

Alouette II Sounder. The lack of significant variations is typical 
of the daytime topside ionosphere outside the plasmapause- The plasma 
scale height indicates the presence of 0* and an if*" flow of 1.5 x 10^ 
ions/cm^/sec can be inferred. 

Figure 15. Plasma Drift Motion Outside the Plasmasphere. The closed curves give 
the trajectories of plasma motion obtained from projecting the 
equatorial drifts of Figure 2 Into the ionosphere. The shaded area 
corresponds regions where the magnetic field lines connect to the 
magnetosheath. Measurements show that the speed of motion varies 
greatly with highest velocities being seen during polar substorms- 

Figure 16, Observations of Plasma Convection. These data represent the observed 
plasma velocities at Chatanika, Alaska (L «* 5*6) on February 11-12, 

1972 at an altitude of 163 km. The dial gives local time and directions 
are related to magnetic north at the center, A substorm is believed to 
have started near 2200 LT and continued through 0400 LT, The overall 
behavior agrees in many respects with the model shown in Figure 15, 
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Figure 17* Reconnection Shock Front Motions* A model of plasma flow along an 
L = 6 field tube has been used to give the downward motion of the 
shock front which attempts to convert the supersonic polar wind to 
subsonic speeds (25]. The shock starts at the magnetic equator and 
slowly moves towards the earth. Because convection of the field 
tube is faster than the shock motion a transition to subsonic flow 
probably does not occur* The flux is the same in either case* 
however. 


DISCUSSION 

Martelli : Could you expand on the assumptions you have used for your theoretical 

model? In particular, what expression for the pressure (scalar? tensor?) have you used 
and what stability criteria for the highly energetic fluxes? 

Banks ; The basic equations are given in a recent paper (Planetary and Space Science 19 
1053 , 1971 ). We assume a scalar pressure and ignore the presence of energetic particles 
Changes in the effective H /0 collision frequency do not greatly affect the H flux 
but do alter the H density profile. 


Hasegawa : What effect does the convective loss of plasmapause particles have on the 

dynamics of filling a flux tube? 

Banks: It is just the convective loss which maintains a low plasma pressure along a 

field tube, thereby inducing' upward flow. 
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Abstract — Vertical profiles of electron density obtain :d in - he vidni ty of the piasmapause using 
the Alouette-II topside sounder have been analyzed to assess the presence of H* flow in the 
topside ionosphere. The observations in the midnight sc: cor show clearly the presence of the 
plasm 
H + flow 

abundance, _ 

to invariant latitudes as iov/ as 43" (L = 2*2) in the regions normally known to be well inside 
the plasma-sphere. The upwards H+ flux is sufficiently large (3 x 10 1 ions sec" 1 ) that 
the plasma -pause cannot be seen in the latitudinal cl :c:ron density contours of the topside 
ionosphere. The cause for this flow remain ; unknown, but it may he a result cf a diurnal 
refilling process. 

1. INTRODUCTION 

During the past several years there has been an increasing appreciation of the dy- 
namical coupling between the ionospheric F s -region and the magnetosphere. In the 
regions exterior to the plasmasphere both the 0 + of the F 2 -layer and the of the topside 
ionosphere are exposed to the general magnetospheric electric field. As a consequence, 
thermal plasma in this region undergoes an £ x B drift which leads to plasma "loss at the 
magnetopause and consequently to the establishment of large scale upward polar wind 
flows of H + . 

For magnetically quiet conditions, the pksmapanse represents the boundary between 
plasma which attempts to co-rotate with the Earth and plasma which is exposed to the 
magnetospheric electric field. (A recent review lias been given by Chappell, 1972.) 
Within the plasmasphere the behavior of H + and Ke + is governed by ion flow along mag- 
netic field lines in response to changes in plasma pressure and chemical composition. 
Such changes occur at low latitudes as a result of F 2 -region ionization and thermal pro- 
cesses and at high latitudes as a consequence cf E x B plasma drift associated with the 
partial penetration of the magnetospheric convection electric field inside the plastnapause 
boundary. The effects of the £ x £ drift arc especially important in the afternoon and 
evening sectors of the plasmasphere where plasma motions across field lines normally 
occur and lead to substantial plasma pressure gradients along the lines of magnetic force. 

The idea of diffusive equilibrium between C : * in the F,- region and the light ions of the 
topside ionosphere has b-eu explored in a number of theoretical studies (Mange, 1960; 
Kockarts and Nicolet, 1963; Bauer, 1966; Hush and Venkateswarcn, 1965). Outside the 
plasmasphere, however, measurements provide evidence for the existence of a much more 
dynamic state with high speed flows occurring at high altitudes acting to replenish regions 
depleted of ionization by £ x B drifts (Bratton a «!., 1971; Hodman, 1969; Banks, 

* On leave of absence from the University of Cttliforain, S.iu Diego, U.S.A, 
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1972). Theoretical studies of these regions (Banks r. .d holier, 1969a, Mmubaflii, 1970; 
Mayr el a/., 1970) have shown the flow to consist principally of moving at speeds near 
or greater than the ion thermal speed. 

The question of light ion flow within the plasmaspherc represents a more difficult 
problem. To a large extent the thermal plasma of the plasmasphcre acts as a fluid con- 
tained in a reservoir with flow along field lines rc:o' ting from changes in pressure or volume. 
A schematic view of the reservoir concept with associated plasma flows is shown in Fig. 1. 
During magnetic storms a large portion of the plasma stored in such magnetic flux tube 
reservoirs is removed by transverse electric fields, either through sweeping away to the 
magnetopause (on the dayside) or through compression and subsequent inward flow to 
tbs ionosphere (on the nightside) [See Chappell (1972) for a more detailed description and 
appropriate references]. Following magnetic storms, ionospheric replenishment of the 
thermal plasma reservoir proceeds via O +/H charge exchange in the dCC-lOQG km region 
and an upward flow of H + into the plasmasphore (Park, 1971; Evans, 1971a, b; Banks 
et al., 1971). Calculations of the thermal plasma behavior for these periods of refilling 
have shown that far out along the magnetic field tubes the H* Sow is Initially supersonic 
and becomes subsonic only after a number of hours have elapsed. In either case there is a 
net outwards H + flow for several days. 



Diffusive Equilibrium: 


Static distribution of plasca with saall invar <1 flov 
at lav eicitudsia lo accwassodate ioo losses In ch* 
F^-rogioti. 



Iav»r d Diffusive Tlavt Plasma ttovej iswari and is last la the F^-rtjzlan.. 

Spaad of f lev cay became supersonic If «xee*?s plasma 
pra**ura Is sufficiently lar#** 



Outward Diffusive Flov: Plasma noves ouhvord sc raualize pressure. 

Sceed of flu*.; riv surer aunlc If ©rersure 

dsflclc Is larne. 


Fis. JL Schematic view of cool plasma behavior in the topside ionosphere. 
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For less is known about the diurnal plasma flow in the absence of refilling effects 
associated with magnetic storms. Incoherent backscatter observations at Millstone Hill 
(L = 3*1) and Arecibo (L — F9) have shown the presence of large (10 s ions cm - - sec _i )0 + 
flows which are believed to be associated, via charge exchange, with H' flow into the 
plasmasphere (Evans, 1971a, b; Ho and Moorcr oft, 1972). Theoretical studies based on 
the equations of continuity and momentum conservation for O* and K + indicate the 
presence of large diurnal variations in H + flow (and consequently G H ‘ flow) resulting from 
changes in temperature and density of the ionosphere (0+) and of the neutral atmosphere 
(O and H) [Nagy and Banks, 1972; Schunk and Walker, 1972; Llayr ?i cil . , 1972; 
Moffett and Murphy, 1973; Geisier and Bowhili, 1965; Hanson and Patterson, 1964]. 
Typical calculated results for the diurnal flux at Millstone Hill are shown in Fig. 2, 
taken from Nagy and Banks (1972). 

The currently available direct observations of topside plasma flow are relatively sparse 
and give an incomplete picture of the important latitude and local time effects. As shown 
in this paper, however, information about flows can be deduced from electron density 
profiles obtained using high altitude radio sounders carried by satellites such as ISIS-I 
and Alouette-IL Such profiles, when interpreted in terms of multi-ion continuity and 
momentum equations appropriate to the topside ionosphere, can provide important data 
on the spatial and temporal extent of thermal plasma flow inside and outside the plasma- 
sphere. 

In the present study, an analysis of Alouette-II electron density profiles has been made. 
The results, discussed in detail in later sections of this paper, reveal the existence of a large 
seals upwards H -1 " flow in the morning sector of the plasmasphere extending from the 
post-sunrise regions into the afternoon sector. In latitude it appears to range from the deep 
polar cap down to A 45° invariant latitude. The cause for this flow is unknown, but 
it could represent either the refilling of flux tubes previously emptied into the night-time 
ivregion or the cross-L drifts of plasma associated with the partial penetration of the con- 
vection electric field into thp plasmasphere. 

The experimental data of this paper are discussed in the next section. The interpretation 



Fig. 2. Computed and observed diurnal H r fluxes at Millstone Hill, Mass accuse nrs 

for 23-24 March 1970. 
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of the data in terms of multi-ion models is described in Section 3, while Section 4 pro- 
vides a general discussion of the results, 

2. MEASUREMENTS 

The Alouette-II satellite includes a sweep frequency ionospheric sounder which permits 
an accurate determination of electron density as a function cf altitude. A complete dis- 
cussion of the equipment and method of analysis needed to convert the apparent range 
frequency records into electron density profiles is contained in papers presented in Free. 
IEEE (June, 1969). 

For the present analysis, a large number of vertical profiles- of electron density were 
provided by the NASA Ames Research Center. These profiles represented data acquired 
principally at Ottawa and Resolute Bay, Canada; vfinkfield, England; and Stanford, 
California. In order to distinguish the effects of hr" flow from elder processes occurring 
in the top side ionosphere it is necessary to examine the electron density profile up to an 
altitude of at least 2000 km and preferably 2500 km. For the present analysis only these 
passes having electron density profiles extending above 2000 km have been used to de- 
termine local time, seasonal, and latitudinal variations in the electron deadly profiles. 

To demonstrate the character of the reduced data, Fig. 3 presents the night-time lati- 
tudinal variation in electron density at fixed altitudes between 500 r.r.d 25C0 k in steps of 
500 k for 7 July 1963. The local time cf this pass extends from 2317 If f at A -= S3 D to 
2329 LT at A — 38° (all latitudes presented are geomagnetic in valiant derived from the 
expression cosVi — XjL where L is the Mcllwain L parameter). The electron density 
contours of this figure show that the variation of electron density with altitude, measured 
in terms of the electron density scale height H c defined as 


1 dn e 
n e dz 




O) 


depends strongly upon invariant latitude. Below the apparent density minimum near 
62°A, the scale heights are relatively large (~ 1300 km) while poleward of this boundary 
the scale heights are smaller (~ 550 km). 



Fig; 3. Contours ojf topside electron density measured by Alouette-II. 

The plasmapauSS is readily apparent In the high altitude densities near 61 °A' while auroral oval 
heating is present between 66 ° and 72°A 




proto.- flow in THE PLASMASEHERE: the morn.no sector 



ftiictron density obtained outside the plvsmaseheee 
Fig. 4. Altitude profiles ° F £ £ rE plas mas?hers (curve 2). 

^SStSSSS « 22 °° kBL ' 

i i_ - u+c- -at fiO° A alo 710 with the simultaneous 
The transition from large to small swl^height - * ’ ;= the topside version of 

arp decrease of electron ^ £’ J and high latitude regions in 

c plasmapause. To show the differences between = ?0 . r and 57.5° on 26 July 

.ore detail, vertical profiles of el ^ tr °” .^ 31 J plotted in terms of geopotential altitudes 
368 are given in Fig. 4 where the dsnsm^ l definition of the scale heignt 

2700 ta a, “” de) *** aris “ fro!I1 a 

issociated with the spatially varying ennn P bot h photoiomzation and plasma 

of 74° A the solar zenith angle becomes le,s man 

heating affect the densities there. important information regarding the dynamic 

Although data suca as t!ls f P r0 ^ p , a5inapa us«, examination of other local 

character of the tops-dc icnospncre ou ^;' . d plasma flows inside the plasma- 

time sectors has revealed the presence of c d u. y r- 1 u i^oortance are upward 
sphere which have not been previously ^ a pat^ah- ^ ^ ^ ^ To 

flows of thermal H" found m tne ^ density profiles obtained in June ana 

illustrate this feature several ^tnadiJ^L ^ ^ ^ local t j rrcs between 07C0 

July 1966 at Stanford anc >r m» ' jcr . s thaa 55^ for alt three passes the maximum 

and 0800 when the solar aenun ang _■ > of 3 July, ^ the maximum 

K at the start of the pass was less tnan 1. tori. 

K in the preceding 7 days was - . 





INVARIANT LATITUDE 


Fig. 7 . Summer morning contours of electron density for 3 July 1966. 
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Fig. 8. Vertical profiles of electron density for 3 July 1966 at three points indi- 

cat ED on Fig, 7, 

These profiles show that 0 + is the dominant to at least 2500 km and, as discussed in th» t« f 
an upward flow of about 3 x 10 3 H + ions cm-* sec" 1 must be present. 


In contrast to the night-time latitudinal profile given in Fig. 3, no discontinuity in 
plasma density or scale height is evident at any of the observed invariant latitudes The 
scale heights are small (300-600 km) and characteristic of O+ alone. There is JittF day to 
day variability present at high altitudes (1500-2500 km) which does not arise from changes 
in the F s -layer densities; i.e. normalization with respect to the 500 km densities tends to 
produce topside densities which are virtually unchanged from day to day Thus it is 
possible to accept the data of those figures as typical for summer morning conditions. 

Typical profiles of electron density as a function of altitude, taken from th“ pass of July 
3 1966, (see Fig. 7) are shown in Fig. 3 for A = 72°, 56° and 50’. When 'the effect of 
gravitational variation is removed, it is found that there is no important change in scale 
height up to 2400 km altitude other than that due to a small increase in plasma tempera- 
ture. As discussed in the next section, such a behavior for the observed conditions can be 
explained only through the presence of an upwards H- flow of magnitude 10* ions cm-- 
sec -1 . 


To show that such a behavior is not directly related to magnetic activity Fig 9 shows 
the latitudinal contours obtained on 9 July 1966 when successive disturbances culminated 
in a K p of 6 at the time of the pass. In spite of this activity, the results of Fi<?. 9 are similar 
to the data of Figs. 5-7 and there is no high altitude plosmapause present to \ - 45° even 
though a jump occurs in the 500 km base density. Note, however, that the plasma scab 
heights are the same on both sides of this change. Suci. a lack of a topside plasma density 
discontinuity provides a distinct contrast to satellite measurements made further awav 
from the Earth nearer the equatorial plane where a plasmapause is normally r . resf , r| » 
(Taylor et al., 1965, 1968; Chappell ct al., 1970, 1971). ' J f 

In order to relate the Alouette-II electron density profiles to the physical processes 
associated with the distribution of ions and electrons in the topside ionosphere it is pee 
essary to introduce the equations describing the behavior of multi-ion plasma’s In the 
following analysis procedure, models of the topside ionosphere arc constructed to match 
the observed profiles. The essential problem is to explain why light ions (H“ or He+) are 
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Fig. 9. Summer morning contours of electron density during disfurbed magnetic 

CONDITIONS, 

Auroral oval heating is present near 70 3 A while variations in the 500 km density are 
reflected in variations at at higher altitudes. 

absent from the summer morning topside ionsphere; i.e. why the typical heavy to light 
ion transition so readily apparent in the night-time data (compare curves I and 2 of Fig. 4) 
is absent in the morning sector (curves 1-3 of Fig. 8). Using the models developed in the 
next section it is found that such a behavior can only be explained if there is an upward 
flow of W and He+. 

3, FLOW EQUATIONS FOR THE TOPSIDE IONOSPHERE 

To analyze the electron density profiles described in the last section we must establish 
a model for the topside ionosphere. Relying upon previous experimental results (Taylor 
etal, 1965, 1968; Hoffman, 1969 and others) and theory we assume that 0 4 is the dominant 
ion in the /^-region and that H 4 is present at higher altitudes, (While tliere is a possibility 
that He 4 may be present, it cannot affect the present analysis since the models must explain 
the absence of all light ions.) 

For the present model we assume that 0+ is the dominant ion at 400 km. Owing to 
diurnal changes in ionization rates, neutral wind momentum transfer, and electrodynamic 
drifts, short ( and long term variations are normally present in the O -1 * density. Thermal 
H 4 * is present in these same regions as a consequence of the accidentally resonant reaction. 

0 + + H ±jZ zz H- + O (2) 

with k f — 2-5 x I0~ n 7^ 2 /j(H) sec -1 and 2*2 x 1 0^ n T} f2 n (O) sec” 1 (Banks and 
Kockarts, 1973). * 

In chemical equilibrium the ratio n(H+)jn(0+) can be obtained as 

*(H+) 9n(H) (T n y* 

n(Q+) 8 n(0) \ Tj ( ) 

so that at 500 km H + is a minor constituent when thermospheric temperatures are greater 
than 800°K. The time constant associated with the equilibrium condition given by (3) is 
short, ranging trom 220 sec in a 750°K thermosphere to T7 sec at 2000'K (Banks and 
Kockarts, 1973). 

Above 700-800 km the condition of chemical equilibrium between O L and becomes 
increasingly difficult to maintain owing to the effects of ion transport arising from plasma 
pressure gradients and gravity. To describe the electron density profiles discussed in the 
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last section, steady state solutions to the equations for the density and flow speed of a 
plasma composed ot and H + are needed. Generally speaking, such a procedure in- 
volves the coupled equations of continuity and momentum conservation for both 0 + 
and (see Banks and Holzer, 1969b). However, since the 0 + flux is much smaller than 
the CH-O limiting flux (~ 2 X IO 9 ions cm -2 sec -1 ), it is possible to ignore the effect the 
° f flow velocity has upon the 0+ density profile (Banks and Holzer, 1969b; Schunk and 
Walker, 1970). In this situation the equation of 0 + momentum conservation can be written 
as 


1 dp(O') 
n(Q+) dr 




(4) 


where s is a coordinate along the magnetic field, p(0 + ) — n(0 + )k T i} a c = n k T , 
= n(0 + ) + n(H + ), gij is the component of gravity parallel to the magnetic field, 
v (0 s H + ) is the 0+-H+ diffusion collision frequency, and r(H+) is the H+ speed along 
the magnetic field. ° 

With regard to H • , fiow effects cannot be ignored and the density and flow velocity 
must be found from the equations 


»(H + ) 


- [n(H+>(H+)] = q(H+) — /f(H+) h(H+) 
dp(H^) 1 dp a 

~~ds h * a dj + w, CH + )g = -m(H + >(H+ 0+) e(H+) 


(5) 

( 6 ) 


where ?(H+) = n(0+MH)*,, /?(H+) = n(0)k b , n(H+) = n(H+)kT u and other quantities 
are similar to those given in Equation (4), 

To obtain the theoretical profiles of ion composition it is necessary to solve Equations 
(4>-(6) subject to given boundary values.. If »(H+) = 0, such solutions describe diffusive 
equilibrium while if c(H- ) 7 = 0, the general flow solutions are obtained describing inward 
and outward motions (Geislsr and Bowhiil, 1965; Geisler, 1967; Banks' and°Holzer 
1969b; Banks et al., 1971 ; Schunk and Walker, i 970). 

For the present study it is necessary to compute ion density profiles consistent with 
the observed electron density and the inferred values of atmospheric parameters. This 
was done in the following way to obtain profiles over the altitude range 500-3000 km: 
the observed electron density at 500 km was taken to be a result of K+ and Cff in chemical 
equilibrium using equation (3). Values of the neutral atomic hydrogen density were 
taken from the satellite results of Mayr and Brinton (1971) for June and July 1366 The 
density of atomic oxygen was inferred from ClRA (1965) models based on the solar 10-7 cm 
radio flux with correction for magnetic activity. 

In addition to this information, the solution to the coupled ion equations requires an 
additional condition relating to the density or flux of H* at the upper boundak As 
shown later, the H+ flux is a relatively poor boundary condition and for this studyaobitions 
to the H+ flow equations were obtained by choosing a specific value of the K : density at 
30C-0 km. Through variation of this upper boundary density it was possible to match the 
ooserved electron density profile with sufficient accuracy to determine the magnitude of 
the flux flowing along the magnetic field lines. 6 

To illustrate this method. Fig. 10 shows the family of H- density profiles appropriate 
to a lOOO^K thermosphere taken from the models of Banks and Kockarts (1973). The 
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Fig. 10. Computed H + and 0 + density profiles for different H j ' fluxes at jOOQ km. 

The O 1 density is largest for the smallest H " profile and vice versa. A lOOCf'X thermosphere 
was used with a 500 km neutral hydrogen density of 1*92 x 10‘ cm" a . The plasma temperature 
(T ( 4- 7’i) at 3000 km was taken as 6500'K. 

electron temperature profile was computed using an electron heat fiux of 5 x 10 y eY 
cm -2 sec -1 with a boundary value T c (500 km) — 3aGG°K and tne usual electron gas 
thermal conductivity. The ion temperature profile was o stained using the ion energy 
balance equations described by Banks ( 1 967) without thermal conduction. For this example 
the magnetic field was assumed to be vertical. 

The individual H+ profiles given by Fig. 10 correspond to different densities at 
3000 lem (10 3 , 5 X 10 3 , 8 X 10-^9 X Iff*, 15 X 10* ions cm~ 3 ). Associated with each 
of these profiles is a computed H+ flux, the magnitude being given directly above each 
curve. Positive values of the flux correspond to outward flows, negative values show inward 


flows. 

The range of O+ density for the different profiles is indicated by the shaded region 
of Fig. 10. When H + is the dominant ion, the 6" density falls off most rapidly with al- 
titude, while when O* is dominant, the slowest decrease is appropriate. This behavior 
results from the electric field coupling between 0“ and M" 7 * and provides an explanation for 
the rapid ion composition changes observed near the plasmapause where Sow becomes 
important (see Banks (1970) for a more detailed discussion of this effect). 

For the purposes of the present study, the most important characteristic of the density 
profiles given ir Fig. 10 is the abrupt change in scale height associated with the transition 
of the major ion species from O* to K + . For tne summer daytime data discussed here, no 
such transition was observable in the electron density p ronlcs up to altitudes of 2000— z500 km. 
Consequently, it is necessary to suppose that the density is low at the upper boundary 
and, as a consequence, that an outward flow is present. 

Although it might be thought that the matching of theory and experiment using this 
method would give only crude estimates of the H+ flux, this is not necessarily the case 
owing to a saturation effect which limits the magnitude oi the flux for a wide range of H ■ 
boundary densities. To show this effect, the value cf the flux has been plotted in Fig. 1 1 
as a function of the 3000 km density using the results, of kig. 10. From tiiis plot we 
find that diffusion equilibrium (he. no H ; flux) occurs when //(il r ) 2DC0 ^3 9 X lG 3 cm 3 . 
For H+ de nsities greater than this value, the flux is inward with a magnitude which increases 
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Fig. I I. Variation of the H flux as a function of the 3003 km H" density. 

The effect of flux saturation occurs quickly an 3 limits the rate at which H ' can leave the 
topside ionosphere. In contrast, there is no limit to the rate a: which H r can Sow into 

the ionosphere. 


■rapidly with increasing H+ density. Such a result indicates that any excess plasma con- 
tained in field tubes above 3000 km can be readily transferred downwards to the upper 
F-region. 

When the H + boundary density is lowered below the diffusive equilibrium value an 
outward H + flux is established. Unlike the inward flow, however, the outward flow rapidly 
saturates such that the H + flux soon reaches a limiting value and further decreases in the 
H + boundary density yield no corresponding change in flux. This behavior illustrates one 
aspect of the limiting flux concept discussed by Hanson and Patterson (1964) and Geisier 
( 1967 ). 


To apply the limiting flux idea to the present study the following line of reasoning has 
been used: the fiux in saturation region depends only slightly upon H f density when 
■the H + density is less than the diffusive equilibrium value. Thus, if wc assume that the H + 
boundary density is equal to the measured electron density at the satellite, we can compute 
a lower limit to the flux. Actually, since there is no change in electron density scale height 
seen at this upper altitude for the present data, the true H + density must be much less than 
the electron density. However, ov/ing to the H + flux saturation effect, the use of the true, 
lower H + density in the computation would not appreciably change the computed value 
of the flux. 

While the above reasoning permits us to compute the H + fiux when the electron density 
profile indicates the absence of appreciable amounts of H + , difficulties arise when the H* 
and 0 + densities become comparable near 2500 km. Reference to Figs. 10 and 1 1 shows 
that the H+ flow is veiy sensitive to small changes in the upper boundary density of H~ 
when H + is the dominant ion above about 1500 km. In this situs, tien small changes in 
the H~ density can lead to large changes in the flux and. consequently, the measured electron 
density profiles cannot be used to predict accurately either the sign or the magnitude of the 
flux. Thus, for example, it is r.ot possible to use the profiles shown in Fig. 3 to compute 
H~ fluxes in the regions inside the nlasmasphere. Although the apparent change in scale 
height shows the presence of H + (or perhaps He 1 ), the uncertainties in temperature, ion 
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composition, and atmospheric densities are sufficient to make any computation highly 

Although chemical and how effects are most important in determining the overall ion 
composition pro i es, tue theimal structure associated with the electron and ion tempera- 

.%L C I n ” 0t ! nCgleaL ' d ; Refereuce t0 Equations (4) and (6) clearly' indicates that values 
temperature^ and weir gradients play a part in determining the final result. For the 
present data direct measurements of temperature were not available from the Alouette-II 
spacecraft. However, the companion satellite, Explorer 31, did carry an electron tern- 

£ r : L ’ rs P rooe wIuch provided (Brace, personal communication) typical latitudinal pro- 
ales of electron temperature for the periods considered here. 

Solutions for the present model require that the electron and ion temperatures be 
Ji%lT Ct iT S of altitude - Because the electron density scale heights were uniformly 
“ ( ^°° km) fo t.^ e l710rnin S passes used here, it was possible to assume that 

, J, aS , C . S f. t0 . a di . ilUSIve equilibrium distribution. Using the experimentally observed 
density distributions it was then possible to solve the equation 

J_ d7; t m(0 + )g H 1 d/I, 

Tp dj ‘ r kT v 1 n t ds ~ 0 (7 > 

for the plasma temperature, T v = T, + T t , by integrating downward from the satellite, 
i he initial value of J* at itie satellite was estimated directly from the slope of the electron 
density scale height at the satellite assuming that 0+ was the dominant ion. Using this 

of ahitude atl ° n ^ ^ ^ USed l ° ° btain the phsira temperature, as a function 

. II to estimate tEle errors in the analysis for the temperature profile. The 

v> ues o T s deouced for the satellite altitude, however, agreed closely with the Explorer-31 
measurements. In addition, changes in composition from O- to H- or He- would give an 
apparent T, higher than true value. Using this method, profiles for T. were deduced 
s owing tnermal gradients of 0-5-1 -CMC knr 1 . Such gradients slightlv modify the H~ 
density profile by making the scale height of H- smaller than would be found for isothermal 
conditions. However, it does not obscure the transition from 0+ to H+ nor charge in any 
important way the H+ limiting flux. 

After thw /alue of T v was determined, computations of T c and T i were made bv assuming 
steady state conditions for the exchange of thermal energy between electrons,' ions, and 
neutrals using rates previously given (see, for example, Banks and Kockarts, 1973). 

4. RESULTS 

Using the mathematical models of the last section, calculations have been made of the 
H flow associated with the typical summer morning Alouette-II pass shown in Fig. 7. 
For tins calculation, the neutral hydrogen density was taken from the results of Mavr" and 
Brmton (1971) The densities of other constituents were determined from standard CIRA 
(19oo) atmospheric models. 

During this summer period the electron density profiles associated with the UoueUe-II 
data give no indication of H-. Thus, to determine die n-fiux.thecomouteJioncomposition 
proules have been made assuming that H+ density at the satellite r.hiiude is equal to the 
observed electron density. As mentioned before, if H~ were really the dominant ion 
near the satellite, then the electron density scale height would be miav times larger 
tian tne values which were observed; the adopted procedure gives an upper limit to The 
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Fig. 12. Computed H + fluxes for 28 June 1966, spanning the range of the normal 

MORNING PLASMAPAUSE. 

These results show that a large upward surge of H + must enter the plasma sphere in the morning 

sector. 


amount of H + actually present and a lower limit to the upward H + flux. Using this method 
the H + vertical flux has been obtained as a function of latitude for the data of Fig. 7 and is 
plotted in Fig. 12. The slight variation of vertical flux with latitude results principally from 
the progressively greater dip angle of the magnetic field. Comparison of Figs. 6 and 13 
shows that changes in the electron density profile with latitude introduce very little variation 
in the flux, and there is no flux discontinuity between the polar wind regions assumed to 
exist above 60°A and the interior of the morning plasmasphere extending down to L — 2-3. 

The results given by Fig. 12 are typical of all morning sector data examined for the 
summers of 1966, 1967 and 1968. However, owing to a lack of atomic hydrogen density 
data no attempt has been made to compute K’“ fluxes for tire later periods. As shown 
elsewhere (Geisler, 1967; Banks and Holzer, 1969b), the limiting H + flux is directly 
proportional to the atomic hydrogen concentration near 6GG-SG9 km. On this basis it 
might be expected that these later fluxes would be substantially smaller than those found 
for 1966. A further difficulty, however, lies in the fact that by decreasing the atomic 
hydrogen concentration one moves the normal, non-flow G+/H+ transition point to higher 
altitudes. Thus, to measure accurately the effects of flow when the hydrogen densities are 
low, one must have data extending to, altitudes higher than the Alcuette-II orbit. 

Examination of the local time behavior of the deduced plasma flow is made difficult 
by the orbital relationships between altitude, local time and season. However, analysis 
of electron density for the summer periods in 1967-1970 indicates that the upward flow 
seems to be present through about 13-14 hr LT. In 1963, however, the neutral hydrogen 
densities were much lower than iu 1966 and the deduction of IF" flow could be made only 
during passes where the electron density profile was obtained up to at least 2500 km altitude. 

5. DISCUSSION 

The H + fluxes described here appear to be a normal feature of the morning sector of 
the plasmasphere. Although magnetic storms are thought to give similar flows, the data 
used in this study show no obvious differences between observations made after many 
magnetically undisturbed days and those obtained during such periods (compare the 
contours of Fig> . 5-7 with Figs. 9). 
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The strength of the flow is sufficiently large to deplete greatly M+ in the topside iono- 
sphere and greatly alter the Q + /H f transition altitude. In this situation it is not possible 
to distinguish the plasmapause at these altitudes, i.e. there is no essential difference be- 
tween the high latitude polar wind and the present flow observed at midlatitudes. Such a 
conclusion may also be supported by a recent study of topside troughs by Tulunay and 
Sayers (1971) where there were essentially no density troughs observed in the mornin° 
sector at Ariel altitudes. However, the Ariel data also reflect changes in composition 
arising from thermal effects. At higher altitudes of Alouette-II, however, the reason for 
such an absence of light ions clearly relates to the lack of H+ or Her inside the plasmasphere 
at iompheric heights; i.e. a flow of light ions is present. Even though the plasmapause is 
not apparent in the topside ionosphere there appears to be no ambiguity in satellite 
measurements made further away from the Earth where the H + density is observed to drop 
to low densities without any indication of a decrease in plasmapause latitude .(Chappell, 
personal communication). This lack of correspondence has not yet been adequately 
explained. 

Incoherent scatter measurements of plasma density and drift velocity can also provide 
information about the H + fluxes described here. At Millstone Hill (L - 3-1), Evans 
(1971a,b) has observed large upwards 0 + fluxes following sunrise and continuing through 
late morning hours. Although a large fraction of this flow is associated with the thermal 
expansion of the 0 L in the topside ionosphere resulting from increases in plasma tem- 
perature, there appears to be a residual component associated with Q+-H charge exchange, 
and subsequent flow. Using Millstone Hill data for 23-24 March 1970 Nagy and 
Banks (1972) made computations of the diurnal H* flux and density arising from changes 
in 0+ density at 500 km. The results indicate an upward flux throughout most of the 
morning. However, the H + density at 3000 km was computed to be about 2 x I0 J cm- 3 , 
about a factor of ten larger than is observed in the same local time sector by Alouette-II. 

Although the existence of the morning I-U Sow appears definite, the causative mecha- 
nism remains to be identified. The flow itself represents the effect of plasma pressure 
gradients along the magnetic field. It is possible that the flow arises as a replenishment 
of ionization previously lost to the F 2 -layer via H — O charge exchange during the summer 
night. As shown by Fig, 1, the plasmasphere undoubtedly acts as a reservoir of ions which 
are able to flow into the A,- region when the 0 + pressure is sufficiently low (Hanson and 
Patterson, 1964). Following sunrise, the /yregion 0+ densities rise rapidly due to 
photoionization and, owing to O+-H charge exchange there is a consequent H+ flow 
upwards into plasmasphere. 

Although the reservoir idea may be partially correct, it ignores the effects of transverse 
electric fields associated with the penetration of the magnetospheric electric field into the 
plasmasphere. oucii electric fields on the morningside oi the plasmasoiiere have a tendency 
to drive plasma away from the Earth (see Sharp ex a 1972) so thata pressure gradient is 
established along the magnetic field driving plasma (H + ) upwards. Such an effect should 
be readily seen in the Tyregicn as a poleward drift of the F-region plasma. This plasma is 
not lost from the Earth, however, since the opposite effect takes piace in the evening 
sector of the plasmasphere wnere plasma is driven inwards (compressed) and downward 
flows of plasma must remove excess ionization from the field tube reservoirs. 

Finally, it is possible that the flow represents a transfer of ions between summer and 
winter conjugate regions (Mayr et al, 1 972). Such a flow could arise from the presence of 
differential pressures between the conjugate F-regions. In addition, the influx of ionization 
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mg. 13. Electron density contours for winter morning conditions. 

There is no evidence for large H + flows from the conjugate region although light ions are 

present below 58' A. 

csn result in substantially greater H + densities in the winter topside ionosphere. In con- 
nection with this possibility, a typical winter morning pass of Alouette-H as seen at Ottawa 
is shown in Fig. 1 3. Although the crowding of electron density contours at low latitudes is 
suggestive of an H + influx, a strong solar zenith angle dependence is also present leading 
to considerable variation in plasma temperature and F 2 -region densities. Examination of 
records of this type has not provided any definite support for the idea of plasma flow 
between conjugate hemispheres. 

From the foregoing discussion, it appears that while the presence of a large scale K + 
flow in the morning sector of the plasmasphere seems well established, its cause remains 
to be determined. 
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Abstract 


Previous statistical studies on the connection between the midlatitude 
ion/electron trough and the plasmaspause have employed a variety of primarily 
non-physical definitions for the trough position and, consequently, have left 
the association of the two regions unclear. In this paper we examine these 
studies and propose a physically motivated definition of the midlatitude trough 
position which, in combination with new information of the filling state of the 
plasmasphere, leads to a clear picture of the connection between the midlatitude 
trough and the plasmasphere. We find that within the plasmasphere the important 
boundary which influences the midlatitude trough is the transition from the 
inner to the outer plasmasphere. This boundary, and the trough, coincide with 
the plasmapause only under special conditions. The inclusion of the inner- 
outer plasmasphere boundary in the general picture of the thermal plqjsma 
distribution in the inner magnetosphere leads to an understanding of the 
relationship between the midlatitude ion/electron trough, the plasmasphere and 
the plasmapause. 



The Physical Connection Between the Plasmasphere and 
the Midlatitude Ionospheric Trough 

C* R. Chappell and P. M* Banks 
Department of Applied Physics and Information Science 
University of California, San Diego 
La Jolla, California 92037 

Introduction 

There has been a great deal of work done to relate the midlatitude 
trough that appears in electron and ion density in the ionospheric F-layer 
to the plasmapause which is observed at higher altitudes in the magnetosphere, 
[Carpenter, 1965, Chappell, 1972], Since the first observations of the mid- 
latitude trough by Muldrew (1965) from Alouette I sounder data and by Sharp 
(1966) from in situ ion measurements, various authors have sought to connect 
the magnetospheric and ionospheric results with a simple physical explanation 
[Thomas and Andrews, 1968; Rycroft and Thomas, 1970; Rycroft and Burnell, 1970; 
Mayr, et al. , 1970]. In these studies limited success has been achieved, at least 

for the nightside region, between local times of 2100 LT and 0600 LT. But if the 

♦ 

underlying assumptions and definitions of these studies are examined closely, 
one finds that the commonality of the ion trough and the plasmapause is far 
from established. And, in fact, one finds that physical processes other than 
those discussed by previous authors may play the dominant role in establishing 
the midlatitude ion trough at a position different from the ionospheric pro- 
jection of the equatorial plasmapause. 

In this paper we first examine the results of the early statistical studies 
of the location of the trough and plasmapause. Following this, we present new 
data on plasmasphere and ion trough characteristics which offer a potential key 
to understanding their association. Finally, we present computed models of the 



topside ionosphere for diffusive equilibrium and ion flow situations to show 
how these models can qualitatively explain the ionosphere and plasmasphere 
observations and provide a link between these two geophysical phenomena* 

Previous Studies of the Midlatitude Trough - Plasmapause Association 

The early studies of the coincidence of the trough and plasmapause loca- 
tions were based on electron and ion density measurements made in the topside 
ionosphere by radio sounders and in situ probes on polar orbiting satellites 
as compared with plasmapause measurements made by the whistler technique [Thomas 
and Andrews, 1968; Ry croft and Thomas, 1970; Ry croft and Burnell, 1970; Muldrew, 
1971; Brace and Theis, 1973]. The regions of the inner magnetosphere involved 
in these comparisons are illustrated in Figure 1. Here we show the inner 
magnetosphere with the ionosphere and plasmasphere as prominent features. The 
specific regions under comparison are illustrated by partial satellite orbital 
tracks with comparison sketches of the expected electron or total ion density 
profiles. (The difference between the electron and light ion density profiles is 
discussed later in this paper.) The dashed curve represents the track of a 
polar orbiting satellite in the topside ionosphere which measures the midlatitude 
trough region as illustrated in the upper sketch. The dotted track covers the 
equatorial plasmapause region with the density profile illustrated in the lower 
sketch. This lower sketch can be derived both from whistler observations and 
direct satellite measurements. Given the existence of the two pertinent sets 
of data it is natural to compare the results and ask in what way do the trough 
and plasmapause compare in position and if they differ, why do they, physically? 

Figure 2 shows a specific example of in situ ion composition data 
measured in the topside ionosphere by the ion mass spectrometer on tigo 4 
(Taylor, Private Communication) . Note the trough locations formed by the 
decreasing H and the increasing 0 (not too pronounced in this case) in the 
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region of about 40 Q to 70° dipole latitude on either side of the magnetic 

•f 

equator. A trough appears very prominently in the H ion density as a decrease 
of about three orders of magnitude through 40° to 70° dipole latitude. A 

» -f- 

transition from H dominance to 0 dominance is also seen at the trough. This 
is a typical effect at the nightside trough location associated with outward 
flow of H + and has been discussed by Binton et^ al. , 1971 (also, in theoretical 
detail, in Banks, 1972). 

The low latitude side of the nightside trough is generally attributed to 
the presence of the plasmasphere while the high latitude side or "cliff’* 

[Rycroft and Thomas, 1970] is attributed to enhanced ionization and higher plasma 
temperatures resulting from auroral precipitation. 

It is important in Figure 2 f to note the very broad range of dipole 
latitudes over which the H + density decreases; i.e., on the low latitude side of 
the trough the initial decrease begins at 30° and extends to 60° dip latitude. 

In L-value this latitude range is approximately from L=1.4 to 4. This spread in 
L value is comparable to the maximum variation in position of the equatorial 
plasmapause during moderate to disturbed magnetic conditions and leaves open a 
broad range of latitude at which the trough position could be defined. It must 
also be noted that the spread of L-shells of the H**" density decrease on the low 
latitude side of the midlatitude trough differs from that seen at the plasmapause 
where the decrease is generally much more pronounced (up to three orders of 
magnitude) and Is confined to a few tenths of an L-shell which is about 2° of 
latitude at L=4. Typical examples of the sharpness of the equatorial plasma- 
pause can be seen in the results of Carpenter, 1966; Taylor et al. , 1970; and 
Chappell, 1972. 

From the foregoing, it is clear that in making a comparative study between 
the ionospheric trough position and the equatorial plasmapause position, one is 
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forced to adopt specific definitions of each, and this is where the problems 
have developed. 

Identification of the plasmapause position outside the ionosphere is not 

3 

difficult, although there has been some discussion over whether the 100 ion/cm 

3 

on the 10 ion/cm density level or the maximum density gradient position should 
be used (see Chappell ert al. , 1970) . The density gradient at the plasmapause is 
generally steep enough and prominent enough in most cases that any of the above 
choices gives essentially the same result within an accuracy of a few tenths of 
an L-shell* 

Such is not the case, however, with the midlatitude trough position. The 

wide range of L-shells covered by^ the trough has led to a number of ad hoc 

definitions. Three definitions employed in earlier trough studies are shown 

in Figure 3. In these examples the troughs were defined in terms of the total 

electron density and therefore represent the latitudinal profile which results 

from the summation of the separate ion densities, i.e., from Figure 2, primarily 
4 1 

H and 0 . Rycroft and Thomas, 1970 (part (a)), chose the minimum electron 

density point in the midlatitude trough to define the trough location. In 

contrast. Miller, 1970 (part (b)), felt that the change in slope at the low 

latitude side of the trough bottom provided a more meaningful definition. 

Brace and Theis, 1973, chose an intermediate point where a straight line 

approximation to the electron density profile on the low latitude side of the 

3 

trough crossed the 1000 electron/cm point. 

Each of these definitions gives a wide variation in the L-shell of the 
individually defined trough positions and, consequently, leads to questionable 
results when comparisons are made with the equatorial plasmapause location. 

Moreover, it appears that the definitions of the trough location are 
consistent only with the general idea that the midlatitude trough must be 
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associated in some way with the density decrease at the equatorial plasma- 
pause. In retrospect, it appears that a more physically significant definition 
may have been overlooked completely. From current theories of ionospheric- 
magnetospheric coupling, one would expect that the physically important point 
in the midlatitude trough profiles is the point where the density begins to 
decrease on the low latitude side of the trough because it is there that 
changes begin to take place in the distribution of ionization in the topside 
ionosphere as a consequence of upward light ion flow. There are certainly 
variations in density at low latitudes associated with neutral winds and 
perpendicular electric fields which may occasionally confuse the initial density 
decrease point. However it is possible in most profiles to identify the point 
where the density begins to decrease, forming the low latitude side of the trough. 
In the data presented from sounder and in situ probe experiments to date this 
point is generally found in the range of L=2 to 3. Statistically, this point 
is approximately two L-shells inside the plasmapause position and this glaring 
difference may have obscured its importance. Whatever the motivating reasons, 
the previous trough definitions are non-physical in nature and consequently have 
led to some confusion in understanding the physical relationship between the 
trough and the plasmapause. As we show here, the use of the point of initial 
density decrease on the low latitude side of the trough can lead to a clear 
picture of the trough and plasmapause association in terms of ionospheric- 
magnetospheric plasma exchange. 

A further clarification should be made about studies of the midlatitude 
trough. This clarification was first pointed out by Taylor and Walsh (1972) 
and concerns the importance of studying the trough signature in light ions as 
opposed to total ion (or electron) density. Taylor and Walsh were able to show 
that the light ion ionospheric trough, as seen in ions, was present at all 
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local times, as is the plasmapause, and therefore provides an improved way 
to study the trough— plasmapause association. This is particularly true in the 
dayside topside ionosphere where the high densities of the dominant photo- 
ionized 0 + ions tend to mask the changes in H + density expected at the ionospheric 
extension of the plasmapause. Although the work of Taylor and Walsh is 
important to studies of the trough it does not alleviate the problem of defining 
the precise trough position since even if one examines H variations (as opposed 
to total ions or n ) , the range of latitudes covered by the decreasing density 
on the low latitude side of the light ion trough is large (see Figure 2). Never- 
theless, if we adopt a trough position using the suggested definition of the low 
latitude point where there is an initial H density decrease, the light ion trough 
position occurs at L-shells which are well inside the statistical plasmapause 
position. An example of this behavior can be seen, for example, in Taylor and 
Walsh, 1972, Figure 5, where the density decreases in light ions are found in the 
range .of 30° to 50° dipole latitude. 

Using the definitions discussed above, previous studies of the trough have 
reached the following conclusions: (1) On the nightside (local times of 2100 

to 0600) the trough position and the plasmapause are close to the same field line 
and move to decreasing L— shells with increasing K.^ [Ry croft and Thomas, 1970; 

Thomas and Andrews, 1968; Burnell and Rycroft, 1970]. (It should be noted that 
a definition of the trough was used which gives locations roughly two L-shells 
poleward of the initial low latitude, light ion density decrease.) (2) Using 
the light ions Instead of the total density, the trough can be seen at all 
local times and is near the expected ionospheric projection of the plasmapause 
(Taylor and Walsh,. 1972). (3) Using total electron density at high altitudes 

(>2500 km where H + is assumed to be the dominant ion) the trough is visible at all 
local times and appears to be circular at about L=4 displaying, in contradiction to 

t-6— 



equatorial plasmasphere measurements, no bulge region (Brace and Theis, 1973). 

From the foregoing, one finds that the trough and plasmapause appear to 
be co-located on the nightside, but, from (3) above, the trough appears not to 
exhibit the prominent bulge region seen in the outer plasmasphere at dusk. 

Although real physical differences may exist between the ionospheric 
projection of the plasmapause and the trough, it is worthwhile to ask whether 
or not the lack of agreement is simply a matter of the ad^ hoc nature of the 
original definitions; i.e., if one uses a more physically meaningful definition 
of the trough, will there be an improvement in our understanding of the 
association between the trough and plasmapause? In the following work, our 
procedure will be to define the trough position as the point of initial decrease 
in light ion density on the low latitude side of the trough. Using this 
definition an examination is made of measurements of the midlatitude trough and 
plasmasphere in the light of current theories of thermal plasma flow. 

Evidence of Disagreement Between the Midlatitude Trough and Plasmapause 

Although we have already pointed out that the statistical studies of the 
trough - plasmapause connection are limited in their successes and challenged 
in their basic definitions, the evidence concerning the disagreement between 
trough and plasmapause position is only indirect. The first study showing this 
disagreement explicity was made by Banks and Doupnik, 1973. (It was also alluded 
to in an earlier study by Br inton jet al. , 1971 in their discussion of a "transi- 
tion region"). Banks and Doupnik examined 'Alouette II sounder data on the 
electron density profiles of the topside ionosphere (up to about 2500 km) near 

dawn. Using models of the topside ionosphere based on the momentum and 

*f 

continuity equation for 0 and H , they were able to interpret topside electron 
density profiles in terms of a large upward flow of H + ions. They found that 
this upward flow was present inside as well as outside the plasmapause and could 
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be found down to invariant latitudes as low as 50° A (L=2.5). At this low 

latitude point there was evidence of a trough in the total electron density. 

But as Taylor and Walsh, (1972), pointed out, this is to be expected because 

of the masking effect of the dominant 0 + in the sunlit ionosphere. Part (b) 

of Figure 4 shows three topside density profiles measured at three different 

latitudes on the pass, A=72°, 56° and 50°. As Banks and Doupnik point out, 

these profiles are consistent with there being a full flow of H ions upward 

out of the topside ionosphere into the plasmasphere. This flow regime is 

generally expected in the plasma trough region outside the plasmapause (Brinton 

et al. , 1971), but it is surprising to see it inside the plasmapause. Only in 

profile (3), which was measured at A=50° (L=2.5) (well inside the expected 

plasmapause position), does the upturn in the density profile near 2000 km occur. 

+ 

This upturn is evidence for a transition to a distribution where H is in 
diffusive equilibrium. 

Xt therefore appears from these data that the transition from a diffusive 
equilibrium distribution to a polar wind— like upward flow distribution can take 
place well inside the plasmapause. Apparently, the transition is not always at 
the plasmapause as has traditionally been thought and this should have an 
important influence on the trough position in the topside ionosphere. There must 
then, be another boundary besides the plasmapause which influences the distribu- 
tion of ionization in the topside ionosphere. 

The Boundary Between the ” Inner and Outer” Plasmasphere 

The explanation for a transition from diffusive equilibrium to full flow 
ion distribution inside the plasmapause can be found by considering the filling 
state of the plasmasphere. In particular, it is necessary to ask if all of the 
flux tubes inside the plasmapause are full of plasma (and therefore in diffusive 
equilibrium) subject only to day-night ebb and flow or if some of the magnetic 



field tubes in the outer portions of the plasmasphere are still in a state of 
filling. If the tubes of the outer plasmasphere have not come to diffjisive 
equilibrium, then there can be large upward flows of ionization present on 
these flux tubes (an idea of the way filling proceeds is given by Banks et 
al. , 1971). 

Data on the degree of plasma filling in the outer plasmasphere has 
recently been presented by Park (1973) . This information is derived from 
whistler data and gives the length of time required for the flux tubes inside 
the plasmasphere to reach a saturation level following a geomagnetic storm. 

Some of Park’s recent data are shown in Figure 5. Here the measured total 
electron content of flux tubes ranging from L=2 out to L=6 are plotted as a 
function of t im e following an isolated magnetic storm which drove the plasma- 
pause to L— shells as low as 2 or 2.5. The storm onset occurred on June 16, 

1965. The curves show the subsequent tube content starting on June 18 and 
continuing through June 25. It is seen that as the tubes fill they eventually 
reach a saturation value of total electron content which increases with increasing 
L as a consequence of increasing flux tube volume. 

An important aspect of these measurements is the number of days required 
for the flux tubes to reach a saturation content. Using correlated bottomside 
sounder data from several stations during this period Park has also found that 
the F 2 *-peak of the ionosphere was able to recover to its average density in a 
period of about three days but it took up to eight days for the plasmasphere 
flux tubes near L=4 to reach saturation, and in this case the plasmapause itself 
is estimated to be outside of L=6. Therefore, we see data that following a 
magnetic storm there is a boundary established earthward (equatorward) of the 
post storm plasmapause inside of which the ionosphere and plasmasphere have 
reached a near diffusive equilibrium distribution and outside of which the flux 
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tubes are still filling with large upflows (Park has estimated the daytime 

8 2 

up flow to be 3 x 10 ions/cm sec at 1000 km with a nighttime downflow of 
8 2 

1.5 x 10 ions/cm sec for this case). As illustrated in Figure 5, this 
boundary between the "inner and outer" plasmasphere expands gradually outward 
(and poleward) to higher L-shells. Note that the boundary is at L=2.5 on 
June 18 (2 days after the storm) at L=2.7 on the 19th, at L=2.9 on the 20th and 
so on. The precise position of the inner-outer plasmasphere transition can be 
quite variable but in general it should depend on: (1) the length of time 

since the previous magnetic storm and (2) the magnitude of the previous storm, 
i.e., the innermost penetration of the plasmapause. . 

Figure 6 shows a conceptual view of the thermal plasma regimes in the 
magnetosphere with the addition of the new boundary of the inner-outer plasma- 
sphere. The innermost or lowest latitude region comprises the inner plasma- 
sphere. Here, the flux tubes are "full" and have reached diffusive 
equilibrium. The outer boundary of this region is initially co-located at the 
storm time plasmapause position and moves outward as the plasmasphere fills and 
the plasmasphere and ionosphere come into diffusive equilibrium (or reach the 
"saturation level" as described by Park). Nevertheless, the inner— outer plasma- 
sphere boundary does not always have to begin at the previous storm-time plasma- 
pause. One can envisage a situation following a large magnetic storm where the 
inner-outer plasmasphere transition zone begins to move outward in response to 
filling during a very quiet period. Then, the onset of a very small magnetic 
storm could move the quiet time plasmapause inwards to lower L-shells but not 
completely into the inner-outer boundary which would continue to fill and move 
outwards. ’ 

The second zone shown in Figure 6 is the outer plasmasphere. Here the 
flux tubes are in a filling state with their plasma distribution in an inter- 
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mediate state between diffusive equilibrium and full flow. In this region 

-f 

there are large upward flows of H each day even though they have significant. 

densities of ions collected near the equatorial plane (See Park, 1973, Figure 3). 

The outer plasmasphere is bounded on the inside by the inner-outer plasmasphere 

transition, and on the outside by the plasmapause. If the global magnetic 

activity remains steady, the outer plasmasphere can gradually fill until the 

inner-outer plasmasphere boundary reaches the plasmapause position. 

Outside the plasmapause the plasma trough is a region of full upward H 

flow similar to the outer plasmasphere with the exception that in this region 

the flux-tubes are emptied approximately once each day in the afternoon-dusk 

sector by the effects of the magneto spheric convection electric field (see 

f 

Review by Chappell, 1972). Although both the outer plasmasphere and plasma- 
trough regions are probably regions of full flow, the plasmapause position 
is still quite evident between them because of the loss of ionization each day 
in the plasmatrough and the steady accumulation of ionization in the -outer 
plasmasphere. 

Finally, the polar cap sector of open magnetic field line is also a full 
upward flow region (see Banks and Holzer, 1969). Here the ion flow is along 
open field lines and is a continual outflow with no accumulation in the equatorial 
regions as occurs in the other zones. 

The Effects of the Inner-Outer Plasmasphere Boundary on the Topside Ionosphere 
The upward plasma flows in the different regions shown in Figure 6 can 
have a significant effect upon ion and electron densities in the topside 

ionosphere. The expected effects are illustrated in Figure 7, taken from Banks 

4 - + 

and Doupnik, 1973. This figure shows the distribution of 0 and H as a 
function of altitude for a variety of flow conditions ranging from a downward 

QJ Q 2 

flow of 2.3 x 10° ions/cin sec at 3000 km to an upward flow of 1.35 x 10 ions/cm 
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sec. The shaded area shows the range of the 0 distribution of the different 

H flows with the upper curve of the shaded area corresponding to tfi'e upward 
+ 8 2 

flow of II of 1.35 x 10 ions/cm sec. The details of the solution of the 
different flow distributions are given in Banks and Doupnik, 1973. Similar 
solutions have also been given by Mayr et al. , 1970 and Brinton et al. , 1971, 
to illustrate the changing topside ionosphere distribution at the plasmapause. 

In Figure 7 the varying flows that are illustrated are a direct result 
of the external plasma pressure on the topside ionosphere. This varying pressure 
is a direct result of the degree to which the flux tubes have been filled and 
how the ionization has distributed itself along the field line during filling 
(see Mayr et al. , 1970 and Banks et al. , 1971, for a discussion of the filling). 
At 2500 km, for example, there can be over an order of magnitude change in the 
density as one goes from a zero upflow (arrow) to an upflow of 1.35 x 10 8 ions/ 
cm sec. Therefore, when a satellite at 2500 km passes from the inner plasma- 
sphere to the outer plasmasphere , it passes from a diffusive equilibrium or low 
flow region to an intermediate filling region which still contains large upflows. 
As a consequence, at this boundary there should be a density drop of an order 
of magnitude or more. 

This effect is schematically illustrated in Figure 8 where the expected 
2500 km electron density is shown for increasing L-shells ranging from the 
inner plasmasphere into the plasmatrough and auroral zone. Within the inner 
plasmasphere the density remains relatively constant (excluding F^-region 
dynamical effects) until the inner— outer plasmasphere boundary is reached. 

Here the electron density drops sharply by over an order of magnitude as a 
result of the transition from H + diffusive equilibrium to a flow regime. 

Owing to uncertainty in the way plasmasphere refilling occurs, we have included 
two dashed curves as alternative ways in which the density fall off might 
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appear with increasing L. Until we can gain further information on the 
actual distribution of ionization along field lines in refilling regions we 
can only speculate on how the plasma pressure in the topside ionosphere should 
vary with L-value in the filling region of the outer plasmasphere. It may be 
possible, however, to use existing satellite data on the measured slope of the 
low latitude side of the trough to gain information about the pressure changes 
in the outer plasmasphere. 

If the outer plasmasphere region is in full upward flow, there should be 

no trace of the plasmapause in the topside ionosphere. If the external pressure 

difference (parallel to the magnetic field lines) present at the plasmapause can 

be transmitted to the topside ionosphere, then there would be a change of flow 

' r 

with a corresponding change in topside density. Finally, the auroral zone and 
polar cap portion of the topside profile is strongly affected by variations 
in F 2 ~layer density ionization caused by electron precipitation. 

The effect of the inner-outer plasmasphere boundary in forming' the low 
latitude edge of the topside ionospheric trough has been verified by at least 
one correlative study of the two regions. Part (a) of Figure 9 shows H 
measurements from an outbound pass of 0G0-5 through the plasmasphere in the after- 
noon sector at 0133 — 0251 UT on August 18, 1973. Part (b) of the figure shows 
an Alouette II pass over the same L— shells about two hours later in universal and 
local time. If one assumes approximate corotation of plasma at these lower L- 
shells, the two passes should be looking at very nearly the same flux tubes. We 
note that in the 0G0-5 data there is a sudden drop in H + density at L=2.6 (giving 
evidence of a previous storm time plasmapause) and another drop at L=6.5 (the 
normal bulge region plasmapause) . The drop at L=2 . 6 probably represents the 
boundary between the inner and outer plasmasphere which appears prominently in 
the Alouette II density data at L=2.5. This correlation verifies the close 
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connection between the midlatitude trough and the inner— outer plastnasphere 
boundary. The density drop of almost two orders of magnitude at L=6.5 in the 
0G0— 5 data, however, is not obvious in the Alouette data. Clearly, in this 
case at least, the density variation which would be characterized as the mid- 
latitude trough in the Alouette II data corresponds to the inner-outer plasma- 
sphere boundary and not the equatorial plasmapause which falls at L=6.5. 

The Relationship Between the Midlatitude Trough and the Plasmapause 

It appears from the above discussion that the initial low latitude density 
decline associated with the midlatitude trough is the transition zone between the 
inner and outer plastnasphere. This boundary is only occasionally found at the 
plasmapause. The most favorable conditions for agreement between the plasmapause 
and the low latitude side of midlatitude trough probably occurs on the nightside when 
downflows from the plasmasphere into the ionosphere are present. In particular, 
at the onset of high magnetic activity when the plasmapause is driven to low L- 
values,the inner-outer boundary of the plasmasphere and the plasmapause are 
established simultaneously at about the same L-shell. A second possible set of 
circumstances leading to co-location of the plasmapause and trough could be after 
a long period of steady magnetic disturbance during which the inner— outer plasma- 
sphere boundary has filled out to the corotation-convection boundary which creates 
the plasmapause. Such a situation would require that the plasmapause position 
remain steady during a complete filling time (about 5 to 8 days). 

What is the significance of the dayside trough? In total ion/electron 
density it probably represents only the decrease in 0 density with increasing 
latitude arising from a decreasing solar zenith angle on the low latitude side 
and the increase in density due to precipitation in the polar cusp on the high 
latitude side. In light ions the dayside trough probably is the signature of 
the inner-outer plasmasphere boundary which can be the same as the plasmapause 
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only for the special filling conditions discussed in the previous paragraph. 

Does one expect a "bulge" in the midlatitude trough which is present in 
the plasmasphere? This does not appear to be the case for several reasons. 
First, the inner-outer plasmasphere boundary is generally found at relatively 
low L-shells (L=2 to 3) where corotation is the dominant influence on the drift 
of flux tubes. The convection electric field should have very little influence 
at these low L-shells and, consequently, the projection of this boundary would 
be nearly circular. Consequently, its signature in the ionosphere as the low 
latitude side of the trough should also be roughly circular. A second reason 
for a circular trough within the topside ionosphere may be associated with the 
fact that if the tubes of the inner plasmasphere did drift to higher L-shells 
in the dusk sector, their volumes would increase and the resulting lowering in 
pressure would induce an upflow of H , thereby eliminating the trough signature 



Summary 

In this paper we have introduced the effect on the topside ionosphere 
of the inner-outer plasmasphere boundary as defined by Park and have suggested 
that it is the dominant factor influencing the low latitude side of the mid- 
latitude ionospheric trough. It appears that the plasmapause and midlatitude 
trough are co-located only for special conditions and that previous statistical 
studies have failed to show this distinction because of somewhat arbitrary 
definitions of the trough position. The agreement between the point of initial 
density decrease on the low latitude side of the ionospheric trough has been 
substantiated in a limited number of cases through the use of 0G0-5 and 
Alouette II correlative data. 
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Figure Captions 


Figure 1 — A sketch of the inner magnetosphere showing the plasmasphere and the 

ionosphere. The dashed line illustrates the path of a polar orbiting satellite in 

the topside ionosphere. A sketch of the expected topside density profile is shown 

in the top right-hand side of the figure. The dotted line illustrates an equatorial 

orbiting satellite pass through the equatorial plasmasphere with the expected density 

profile shown in the lower right-hand corner of the figure. The shaded region rep- 
resents the plasmasphere. 

Figure 2 — An example of ion density measurements in the topside ionosphere as 

measured by the mass spectrometer on the Ogo 4 polar orbiting satellite. The trough 
position is evident on each side of the equator at dipole latitudes of 40° to 70°. 

Figure 3 — Three different definitions of the midlatitude ion trough position. 

Part (a) is taken from Rycroft and Thomas (1970) and uses the minimum point as the 
trough definition. Part (b) is from the work of Miller (1970) and uses the change in 
slope as the definition of the trough. Part (c) from Brace and Theis, 1973, uses 
the 1000 electron/cm^ level as the trough position definition. 

Figure 4 — Alouette II data taken from Banks and Doupnik (1973). Part (a) shows 

the density versus invariant latitude for several altitudes on a dawn (0821 local 
time) pass. Part (b) shows three electron density profiles in the topside ionosphere 
measured at three different invariant latitudes along the satellite pass, 72°, 56* 
and 50°. 

Figure 5 - Whistler measurements from Park, 1973, on the refilling of the 

plasmasphere following the magnetic storm of June 16, 1965. This figure shows the 
total electron tube content as a function of L for successive days beginning June 18, 

2 days following the storm. The tubes are seen to fill up to a saturation level shown 
by the heavy solid line. 
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Figure 6 - A sketch showing the different zones of thermal plasma density in 

the magnetosphere. In this sketch the inner plasmasphere has been added to the 
previously described regions of the outer plasmasphere, the plasmatrough, and the 
polar wind. 

+ + 

Figure 7 - A plot of the distribution of H and 0 ions In the topside 

ionosphere taken from Banks and Doupnik, 1973. The different distributions represent 

changing flux conditions at the 3000 km boundary level. The H fluxes at 

8 2 8 2 
3000 km vary from 1,35 x 10 ions/cm sec upward to 2.3 x 10 ions/cm sec down- 
ward. The shaded region represents the variation in the 0 profile for varying H 
fluxes with the top curve of the shaded region representing an upward H + flux of 
1.35 x 10° ions/cni sec. The diffusive equilibrium or zero flux H profile is shown 
by the arrow. 

Figure 8 - A sketch of the expected ion density variations in the topside 

ionosphere. The density drops significantly at the boundary between the inner and 
outer plasmasphere. The exact nature of the drop with increasing L is unknown 
because of the unknown filling distribution of ionization along the magnetic field 
in the outer plasmasphere and plasmatrough. The high latitude increase in ionization 
is due to auroral precipitation. 

Figure 9 - OGO-5 and Alouette II correlative data taken on August 18, 1973 

showing the inner-outer plasmasphere boundary signature in both the plasmasphere 
and topside ionosphere. In this refilling situation this boundary is seen to occur 
at about L=2.6 well inside the plasmapause boundary at 1=6.5. 
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